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We consider an ensemble of spherical colloidal particles immersed in a near-critical solvent such as a binary 
liquid mixture close to its critical demixing point. The emerging long-ranged fluctuations of the corresponding 
order parameter of the solvent drive the divergence of the correlation length. Spatial confinements of these 
critical fluctuations by colloidal solute particles, acting as cavities in the fluctuating medium, restrict and 
modify the fluctuation spectrum in a way which depends on their relative configuration. This results in effec- 
tive, so-called critical Casimir forces (CCFs) acting on the confining surfaces. Using the available knowledge 
about CCFs we study the structure and stability of such colloidal suspensions by employing an approach in 
terms of effective, one-component colloidal systems. Applying the approximation of pairwise additive CCFs 
we calculate the radial distribution function of the colloids, which is experimentally accessible. We analyze 
colloidal aggregation due to CCFs and thus allude to previous experimental studies which are still under 
debate. 

PACS numbers: 61.20.Gy, 64.60.fd, 64.70.pv,64.75.Xc, 82.70.Dd 



I. INTRODUCTION 

In colloidal suspensions with near-critical solvents, col- 
loidal particles acting as cavities impose boundary con- 
ditions at the colloid surfaces for the fluctuating order- 
parameter (OP) of the solvent and perturb the OP 
field on the length scale of the bulk correlation length 
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OJ ~ \t\ v ', which diverges 

(s) 

upon approaching the bulk critical temperature T c of 
the solvent^ v is a standard critical bulk exponent. 
These modifications of the OP and the restrictions of 
its fluctuation spectrum result in an effective force act- 
ing between the colloids, known as the so-called critical 
Casimir force (CCF)»i For equal boundary conditions, 
i.e., equal preference of all colloids for one of the two 
species of the binary solvent, the CCF is attractive and 
can lead to reversible aggregation phenomena. Aggrega- 
tion phenomena are of considerable interest due to their 
relevance for the stability of colloidal suspensions. 

In our preceding paper— we have investigated the phase 
behavior of colloidal suspensions with near-critical sol- 
vents and have explored the possibility of a phase sep- 
aration into two thermodynamic phases, one being rich 
and the other being poor in colloidal particles. If the ef- 
fective attraction potential between colloidal particles is 
sufficiently strong, a condensation of colloids into a stable 
equilibrium liquid phase may be preempted (on the time 
scale of observation) by the formation of non- equilibrium 
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aggregates in which the colloidal particle stick together. 
In general, these aggregates may grow and shrink and 
their structure depends on the packing fraction of the 
colloidal particles in the aggregates and on the strength 
of the attraction among the colloidal particles. In gen- 
eral the aggregates may consist of crystalline structures, 
gels, glasses, but also loose fractal aggregates may occur. 
In the present paper we shall address the possibility of 
aggregation of colloids immersed in near-critical binary 
liquid mixtures serving as their solvent. 

Beysens and Esteve^ have reported aggregation phe- 
nomena for silica spheres immersed in a binary liquid 
mixture of water and lutidine which exhibits a miscibil- 
ity gap with a lower critical point Using light scat- 
tering these authors found that upon approaching the 
bulk coexistence region of demixing from the one-phase 
region of the binary liquid mixture at constant com- 
position of the solvent, the silica spheres coagulate re- 
versibly at rather well defined temperatures. Moving 
back the thermodynamic state deeply into the one-phase 
region the coagulated aggregates dissolve again. These 
observations, i.e., the formation of colloidal aggregates 
which sediment, have been confirmed later for silica and 
polystyrene particles immersed in the same solvent^— 
and have been reported also for other systems jiiilir— for 
reviews see Ref . Il7l . The experimental findings have been 
interpreted as (the- precursors of) a phase transition in 
the ternary mixtureii~— or as a (non-equilibrium) co- 
agulation, which can be caused by various mechanisms. 
Some aspects of these experimental observations will be 
discussed in more detail in Sect. [VTl 

Generically, colloidal particles have a preference for one 
of the two components of the binary solvent. At the sur- 
face of the colloid this preference gives rise to a surface 
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field conjugate to the order parameter at the surface and 
thus leads to an adsorption layer rich in this component. 
Upon approaching the critical point along the two-phase 
coexistence curve, at the wetting transition the adsorp- 
tion layers on the surface of the spherical colloids be- 
come thick (but stay finite) so that they can snap if the 
colloids get sufficiently close.— This leads to attractive 
interactions 1 ^ and can give rise to coagulation. 

Close to the critical point CCFs provide an additional 
mechanism for mutual attraction between the colloids. 
Experimentally it has been observed that (strong) coagu- 
lation occurs for compositions of the solvent for which the 
concentration c a of the adsorbed component preferred by 
the colloids is comparable but less than its value at the 
critical point, i.e., when c a < c a s ,c- This observation is 
in line with both explanations mentioned above because 
the wetting transition is indeed located in that region of 
the bulk phase diagram of the solvent and because the 
CCFs are stronger therein than for c a > Ca; c - Several 
theoretical approaches have been used in order to study 
this aggregation phenomenon for colloids immersed in a 
solvent which is close to its phase separation^ - — Most 
of them, however, stay on a qualitative level of descrip- 
tion, pointing out the possibility of coagulation. 

Here we present a quantitative analysis of the stabil- 
ity of monodisperse spherical colloidal particles immersed 
in a critical binary solvent. Our analysis is based on 
controlled approximations which we have worked out in 
Ref. 0. Following Ref. we adopt a description in terms 
of an effective one-component system. Within this ap- 
proach the presence of the solvent is taken into account 
via an effective pair potential of the colloids. For the 
regular background interaction potential we consider a 
generic, softly repulsive contribution. The CCFs depend 
on the configuration of all colloidal particles and thus 
are in general non-additive. We focus on dilute suspen- 
sions with low colloidal number density p, for which the 
approximation of pairwise additive CCFs is valid. The 
critical Casimir pair potential between two spherical par- 
ticles immersed in a binary liquid mixture is described in 
terms of a universal scaling function corresponding to the 
Ising universality class with symmetry-breaking bound- 
ary conditions! 19 ' 2 5 ~ — We apply the Derjaguin approx- 
imation, within which the pair potential between two 
spheres is expressed in terms of the scaling function of 
the CCF in a slab^r— The Derjaguin approximation 
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In Eq. (QJ the background interaction potential between 

(s) 

the colloids, which is present also away from T c , consists 
of a hard core repulsion for ccnter-to-ccntcr distances 



can be applied, if the surface-to-surface distances D be- 
tween the objects on which the CCF acts, i.e., D < £, 
are small compared with the radius R of the colloids. 

In d = 3 and for a vanishing bulk ordering field hb = 0, 
the corresponding scaling function is known from Monte 
Carlo simulations j^ 3 . - — from field-theoretical studies^ 
and from the extended dc Gcnnes-Fisher local-functional 
method. 38 ' 39 Corresponding experimental dat a 32 ' 40 " — 
arc also available. The conjugate ordering field hb of 
the order parameter of the solvent corresponds to the 
chemical potential difference of the two species forming 
the binary liquid mixture. The reversible aggregation 
phenomenon of colloids is observed for thermodynamic 
states of the solvent which correspond to a nonzero hb- 
In order to handle the dependence of the CCF on this or- 
dering field within the Derjaguin approximation, we use 
still another approximation for the film scaling function 
of the CCF in d = 3 (see, c.f., Eq. ©). The scaling func- 
tions resulting from this approximation are comparable 
with the ones obtained using the extended de Gennes- 
Fishcr local-functional method^ 3 - and are in qualitative 
agreement with corresponding curves provided in Ref. [l6| 
which are based on a density functional approach. 

Our presentation is organized as follows. In Sect. |TT] 
we shortly summarize the effective one-component de- 
scription which we use for our analysis. For more details 
and a comprehensive discussion of the validity and of 
the limitations of this description we refer to the preced- 
ing paper.— In Sect lllll we discuss the possible functional 
forms of the effective interaction potential employed in 
our study. Section IIVI is concerned with aggregation 
phenomena in colloidal suspensions and in Sect. [V] re- 
sults for the radial distribution function of the colloids 
are presented. In Sect. IVII we compare available experi- 
mental data with our results. We summarize our study 

in Sect, eld 

II. EFFECTIVE INTERACTIONS 

We aim for a theoretical description of the behav- 
ior of monodisperse colloids immersed in a near-critical 
solvent. To this end we consider a monodisperse, ef- 
fective one-component system of colloids with radius R 
for which the presence of the solvent enters via the fol- 
lowing effective pair potential V (r = D + 2R) / (ksT) = 
U(r = D + 2R): 



(1) 

I 

r < 2R and a soft, repulsive tail U rep for D = r — 2R > 0. 
The latter prevents aggregation due to omnipresent effec- 
tively attractive dispersion forces. In the case of electro- 



U(r) 



(d=3) 



D < 

Aexp(-nD) + (1/A)6»( d = 3 )(6, A, £), D > 0. 
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static repulsion, the range k _1 of the repulsion is associ- 
ated with the Dcbye screening length and its strength A 
depends on the colloidal surface charge density . 44 > 45 We 
do not consider an additional interaction which would 
account for effectively attractive dispersion forces. Ef- 
fectively, dispersion forces can be switched off by us- 
ing index-matched colloidal suspensions. As will be dis- 
cussed in the next section, the presence of attractive dis- 
persion forces does not change the conclusions of the 
present study. The critical Casimir potential Uc d ~^ 
has a scaling form with 9^ d= ^ as its universal scaling 
function. 9^ depends on the spatial dimension d and 
on the scaling variables 9 = sgn (t) D /£, A = D/R, 
and E = sgn (hb) £ / . The bulk correlation lengths 

£(*£0) =£\t\- v and = d H) \h b r /m govern 
the exponential decay of the bulk solvent OP correlation 
function for t — > ± at hb = and for hb — > at t — 0, re- 
spectively; one has t = (Ti s) -T)/r c (s) for a lower critical 
point. The amplitudes £ ' are non-universal but related 
via universal amplitude ratios; v , /?, and S are standard 
bulk critical exponents 4£ 

For 0( d=3 ) we use the Derjaguin approximation. For 
the configuration of two spheres it renders^— 



/oo 
( 



x~ 2 - x~ 3 ) i?[ d ~ 3 ^ (xQ, E) dx; 

(2) 

i?jj d ^ (y, E) is the universal scaling function of the CCF 
fc,\\ P er ksTc an d per area A for a slab of the thickness 
L: f CA /A = k B T { c s) L- d ^ ] (y = agn(t)L/£,Yl). From 
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Eq. ([2]) it follows that the scaling function 
depend on A, which therefore enters into Eq. ([I]) only as 
a prefactor. In the cases in which E = 0, i.e., for a van- 
ishing bulk field hb we use (y, E = 0) as determined 
from Monte Carlo simulations. 33-36 



The scaling function for nonzero E is obtained from 
the following approximation (see Eq. (24) in Ref. [2| and 
the subsequent discussion there): 



(y,E = o) 



j(<*=4) 



tf[, d = 4) (i,,S = 0) 



(3) 



Within this approximation, for hb — > the scaling func- 
tion reduces exactly to $y (E = 0) (for all d) and for 
fixed values of y its functional form is the one obtained 
from mean-field theory (MFT, d — 4). The neces- 
sary input for this approximation is the mean-field scal- 
ing function of the CCF for the film geometry, which 
follows from the field-theoretical approach within the 
framework of the Landau-Ginzburg theory. For col- 
loids strongly preferring one and the same of the two 
species of the binary liquid mixture the so-called strong 
adsorption limit applies. Accordingly, the appropriate 
boundary conditions are the ones breaking the symme- 
try of the OP at the colloid surfaces. We determine the 
CCF from the local OP profile (j) M FT within MFT via 
the so-called stress-tensor 7"4£ The CCF per area in a 
slab, which is confined along the z-direction, is given by 
the (z,z) component of the stress tensor, f^, FT ^/A = 

k B T^ s) (T z ,z\4>mft] - T z ,z [4>b,MFT]) (see Eqs. (20) and 
(25) in Ref.0). 

Our results will be expressed in terms of the parame- 
ters 



s = kR, 
and x 



a = A/s, 
kD = nr 



C = sgn (t) «;£, 



2s. 



(4) 



Within this parametrization the functional form of U 
is determined by a, £, and E, while s tunes the overall 
strength of the potential without affecting its shape: 



U (r = D + 2R) = 



\ oo 
s 



c x < 

{aexp(-x) + (l/x)<= 3) (x/C,E)} x > 0. 



(5) 



The parameter £ can be varied experimentally by chang- 
ing the temparature (via £) or by adding salt to the so- 
lution (via k); it expresses the competition between the 
length scales of electrostatic repulsion and of the attrac- 
tive CCF. The strength of repulsion is measured by a; it is 
usually kept constant. E is determined by the thermody- 
namic state of the solvent. The ranges of the values of the 
parameters entering into the effective potential (Eq. (JSJ) 
and of the scaling variable E, which correspond to pos- 
sible experimental realizations, are discussed in detail in 
Ref. 0. The divergence ~ a; -1 in U (D = x/k — > 0) (see 
Eq. (|SJ|) is unphysical^ because the scaling form used for 



Uc applies only for distances large compared with mi- 
croscopic length scales such as £^ . In order to circumvent 
this unphysical behavior we shall apply a linear extrapo- 
lation scheme for U {D —¥ 0) whenever this is necessary. 



III. SHAPES OF THE EFFECTIVE POTENTIAL 

By varying the thermodynamic state of the solvent 
the range and the strength of the attractive CCFs are 
changed. This variable attraction combined with the re- 
pulsive background contribution results in a broad va- 
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x = kD 

FIG. 1. Exemplary shapes of the effective pair potential given 
by Eq. The dashed line represents the purely repulsive 
part of the potential for a — 20. The potentials considered 
here are physically appropriate for x = kD > 1. At smaller 
values of x (grayish area) additional (effective) forces can oc- 
cur. In particular, the occurrence of a first minimum at con- 
tact, i.e., at D = 0, depends on microscopic details, which are 
not taken into account by Eq. ((5]). 



riety of possible shapes of the effective potential given 
by Eq. (|5|). Examples of various shapes of this effective 
potential, as obtained for different values of the parame- 
ters a, C and E, are shown in Fig. [l] where we plot the 
reduced potential u = U/s versus the reduced distance 
x = kD. Depending on the values of the parameters, 
the CCFs compensate the repulsion for all values of D or 
only within certain regimes. Compensation for all values 
of D occurs if the strength a of the repulsion is not too 
strong. This is illustrated by the blue curve in Fig.[T]cor- 
responding to the parameter values a = 10, ( = 2.5, and 
E = —4.2. We recall that negative values of E correspond 
to solvent compositions which are poor in the species pre- 
ferred by the colloids. For larger values of a, a secondary 
attractive minimum (i.e., u < 0) can occur (in addition 
to the primary global minimum for D — > 0, which will 
be discussed below), while for small distances D the po- 
tential maintains a repulsive part, as for the orange and 
green curves in Fig.[TJ Since s = kR can vary between ca. 
10 and 100 a shallow minimum in u can correspond to a 
deep potential well in U. For suitably chosen parameters, 
such as for the green curve (a = 20, £ = 2, E = —3.7) in 
Fig. [TJ the potential may be repulsive at large distances 
D — > oo and may have an attractive secondary minimum. 



Thus the CCFs provide a mechanism for realizing short- 
ranged attractive and long-ranged repulsive interaction 
pa ir potentials as studied for example in Refs. H3 and 
|49| in related contexts. These studies have shown that 
such potentials give rise to complex phase behaviors of 
colloidal suspensions. 

The issue of the contact value U (D = 0) and the pos- 
sibility of the occurrence of a primary minimum at small 
values of D is subtle and, as pointed out in Sec.|Hl cannot 
be discussed in general in terms of the effective potential 
U given by Eq. ([5|). However, U captures the occurrence 
of the potential barrier for D ~ n^ 1 and its reduction due 
to the attractive CCFs which is sufficient for the focus 
of the present study, i.e., concerning the influence of the 
universal CCFs. The short-ranged interactions acting at 
distances D <C (with £ 3> £0) depend sensitively 

on materials properties. But for the following reason they 
do not affect our consideration of the stability of col- 
loidal suspensions. In the case of large a, the repulsion 
is strong enough such that the CCFs cannot significantly 
reduce the associated potential barrier. Therefore within 
experimental observation times only the pair-interaction 
potential U considered in Eq. ([5]) is relevant for determin- 
ing the phase behavior and the coagulation process. The 
influence of attractive short-ranged forces could be ob- 
served only for very long waiting-times because the time 
scale for passing the barrier is ~ exp (U) (see Subsect. lIVI 
below). On the other hand in the case of small a the 
CCFs fully compensate the repulsion. The estimates for 
the values of the parameters for which the repulsive bar- 
rier disappears and thus rapid coagulation sets in can be 
obtained from the effective potential U given in Eq. (|5|). 
However, in order to be able to determine the detailed 
process of coagulation a truly microscopic description is 
needed. This is beyond the scope of the present study 
devoted to universal features. As long as no additional 
microscopic forces compete with the considered effective 
potential, Eq. (|5|) provides a description which is suffi- 
cient for the present purpose. 

IV. STABILITY 

The stability of colloidal suspensions and the aggrega- 
tion of colloids are related to kinetic processes (see Ref.l45l 
and references therein). They are based on the diffusion 
of single particles in the presence of other particles of the 
same kind, interacting with them via interaction poten- 
tials which contain both attractive and repulsive contri- 
butions. For interacting particles, which irreversibly stick 
together once their surfaces touch each other, Fuchs in- 
troduced, using Smulochowski's theory, the stability ratio 

W = 2R fggjpO> d , ( 6 ) 

J2R r 

For hard spheres W = 1, while for W > 1 (W < 1) the 
repulsive (attractive) part of the interaction potential U 
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FIG. 2. Contour lines of constant stability ratio W — 1 
(Eq. (|6])) for the effective pair potential given in Eq. ((5} as 
function of a and E for various (fixed) values of f. The star 
and the dagger superscripts correspond to s — 10 and s = 50, 
respectively. For values W < 1 (the region below the corre- 
sponding contour line) the colloids coagulate rapidly. For a bi- 
nary solvent with compositions which are smaller, i.e., slightly 
poorer in the component preferred by the colloids, than the 
critical composition, i.e., for E < 0, this region (W < 1) 
broadens towards larger values of the parameter a of the re- 
pulsion strength. The logarithmic plot in the inset shows, 
that for rather large negative values of E the contour lines 
decay according to a power law as function of Az = z m — z. 
The slope of the dashed line is 0.84; z m = —2.6 approximates 
the position of the maxima of the contour lines. 

dominates. If there is a potential barrier, i.e., U (r) 1 
for a certain range of distances r > 2R, which leads to 
W > 1, on intermediate time scales the suspension will 
equilibrate into a (meta) stable homogeneous state. Only 
on very large time scales cluster formation sets in slowly. 
It turns out that W is proportional to the ratio between 
the characteristic times for diffusion of a single particle, 
tdiff, and for the formation of a pair of particles, t pa ir'^ 

tdiff/tpair = 3r]/W, (7) 

with the packing fraction r\ = (Att/3) R 3 p; p is the num- 
ber density of the colloidal particles. 

We have calculated the stability ratio W = 
W (a, s, C, E) according to Eq. © by using Eq. ([5]) within 
the one-phase region. Since the stability ratio W depends 
exponentially on the potential, in the calculation of W 
we can afford to omit the use of a cut-off at small dis- 
tances D. The main contribution to W stems from the 



maximum of the potential. The use of a cut-off would 
only slightly change the value of W but it would require 
an additional specification which is not necessary for this 
analysis. 

In Fig. [5] the contour lines W — 1 are plotted in the 
space spanned by the parameters a and z = E/£ for 
various values of £. Changing S for fixed £ and s cor- 
responds to moving along an isotherm in the thermody- 
namic space of the system. In order to have the contour 
lines on the same scale, it is useful to introduce the re- 
duced variable E/£ as abscissa. For W = 1 the time 
scales of diffusion and pairing of two particles are com- 
parable (Eq. 0). For a given value of (, in the region 
W > 1 above the corresponding contour line the suspen- 
sion initially equilibrates into a (metastable) fluid state. 
If the potential has its global minimum at D = 0, e.g., 
due to van-der-Waals forces, aggregates may form even- 
tually at long times; these time scales depend on the value 
of W. On the other hand, for W < 1 rapid coagulation 
sets in. 

The overall characteristics of the shape of contour lines 
W = 1 remains the same for different (but fixed) values of 
( and s (see Fig. [2]). Upon decreasing S towards slightly 
negative values (i.e., for a binary solvent upon decreasing 
the amount of the component preferred by the colloids 
relative to the critical composition) the colloidal stability 
remarkably decreases, i.e., the contour line W = 1 is 
shifted to larger values of the parameter a of the repulsion 
strength. For even more negative values of E, the system 
is driven far away from criticality, so that the attraction 
due to the CCFs weakens and stability is favored again. 
The position z max of the maximum of the corresponding 
line a = a (z = E/£) (given by W (a, s, (, E = zQ = 1 for 
fixed s and fixed £) does hardly vary with £. However the 
amplitude of a (z) and in particular the maximum value 
o-max = a{z m ax) exhibit a significantly nonmonotonous 
dependence on ( ~ \t\ ~ v . The value of a max is largest 
for ( ~ 2. This behavior reflects the fact that for (+, +) 
boundary conditions (i.e., equal surface fields on both 
surfaces) as considered here the CCFs are not strongest 

(s) 

at Tc ' and hf, = 0, but slightly within the one-phase 
region and for small negative values of hb- 

The contour lines corresponding to a fixed value of s 
and various values of £ merge to a single curve for large 
negative and large positive values of E/£. This behav- 
ior can be explained by the fact that for large values 
of |E| the scaling functions (© = D/£,E) plotted 

as a function of OE for various values of E fall onto 
each other4£ Interestingly, within the plotted range of 
large negative values of z = E/£ the (common) decay 
of the contour lines can be described by a power law 
~ (Zm ^ z ) P = Az _p (see the double logarithmic plot 
in the inset of Fig.[2l where z m = —2.6 approximates the 
various z rnax (£) andp ~ 0.84). In the limit of strong bulk 
fields, i.e., \hb\ — > oo, the CCFs are vanishingly small. In 
this limit, W = 1 for HS without a soft repulsive contri- 
bution and therefore (a (\z\ — > oo)) w=1 — > 0. The level- 
ing off of the contour lines W = 1 close to E = is a 
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consequence of the choice of the abscissa ~ hb"'^ ■ As 
a function of hb the contour lines exhibit no leveling off 
around 2 = 0. 

The potential under consideration (within the Der- 
jaguin approximation) is proportional to the parameter s. 
Thus with increasing (decreasing) values of s the region 
W < 1 of aggregation shrinks (widens), because the re- 
pulsive barrier in U becomes larger (smaller). However, 
since the repulsive barrier contributes exponentially to 
W, the boundary of the region W < 1 is mostly deter- 
mined by the sheer occurrence of the repulsive barrier 
(which, however, is independent of s) and therefore the 
dependence of the contour line W = 1 on s is weak. This 
is seen in Fig. [2] where for £ = 2 we compare the contour 
lines W = 1 for s = 10 and s = 50 (the blue and gray 
curves, respectively). Note, that this weak dependence 
of the stability region on the parameter s = kR is in 
contrast to the strong dependence on s of the two-phase 
region of the effective colloidal system, as discussed in 
Ref. H (see especially Fig. 3 therein). For the two-phase 
region not only the existence but also the actual strength 
of the attraction is important, and therefore it is sensitive 
to the value of s, because U ~ s. 

The contour lines in Fig. [2] of constant stability ra- 
tio W = 1 can be compared with the contour lines in 
Fig. 3 in Ref. of constant critical temperature tJ 6 ^ 
(or £c ^) delimitating the " liquid "- "gas" coexistence of 
the effective one-component colloidal system. In both fig- 
ures the contour lines are plotted in the plane spanned by 
the parameter a, measuring the strength of the repulsion 
of the effective pair potential, and the scaling variable S 
(or E rescaled by the fixed value of f). The gross fea- 
tures of these contour lines are similar but the values of 
a are rather different. (Note the different scales for a 
in these two figures.) For example, for the same value 
of s = kR = 10, at the temperature corresponding to 
( = 2 strong coagulation (i.e., W < 1) occurs for a < 13, 
whereas at the comparable temperature corresponding to 
C = 3 the colloidal suspension starts to phase separate al- 
ready for a much stronger repulsion, i.e., for a < 300 (see 
the blue curves in Fig. 3 in Ref. 0). We note, however, 
that phase separation occurs at relatively high packing 
fractions r\ > 0.06 of the colloidal particles (see also Sub- 
sect. Ill B in Ref. 0) whereas aggregation occurs also for 
very small packing fractions of the colloidal particles. 

The experimentally relevant case of the variation of 
the stability ratio W with temperature at constant OP 
is shown in a semi-logarithmic plot in Fig. [3] in terms of 
£ ~ \t\ and the parameter too which is related to the 
OP; for t > it is given by (see Eq. (27) in Ref. [3) 

mo = sgn(4>)(C+) 1/w |B/0| 1/0 (8) 

and C + = <o ■ 

(j) and hb are related by the equation of state, hb = 

Vsgn{(t>) \4>\ S F(t\B/<j>^ ,P \ where F (x^j = 1+X (in 
lowest order in X) is a universal (bulk) scaling func- 




2 4 6 

c 

FIG. 3. Semi-logarithmic plot of the stability ratio W 
(Eq. ©) for the typical experimental situation of varying 
the reduced temperature \t\ ~ a ^ g xec j composition 

Ua - <k£\ ~ mo" 1 (Eq. ©). The steep decrease of In (W) 
towards negative values upon increasing f, i.e., for t — > 0, 
suggests that coagulation sets in within a narrow tempera- 
ture interval. The increase of W for large values of ( occurs 
for strongly repulsive contributions to the effective potential 
(a = 10, Eq. (0); s — kR, where At -1 is the range of the 
repulsive "background" potential and R is the radius of the 
colloids. 

tion (see Appendix A in Ref. 0). B and V arc non- 
universal amplitudes which depend on the definition of 
4> such that on the coexistence curve the bulk OP fol- 
lows 4>b (t — > 0~~,hb = 0) = B \tf and that on the criti- 
cal isotherm 4> b (t = 0, hb — > 0) = sgn (h b ) \hbjT>^^ . For 
example, in binary liquid mixtures the order parame- 
ter (j) is proportional to the deviation of the concentra- 

(s) 

tion c a of the component a from its critical value Ca,c, 
4> = A(c a — Ca S /c) (note that B is proportional to A), 
which can be easily controlled by changing the mass or 
the volume fraction of one of the components of the mix- 
ture. 

The steep decrease of W upon increasing £ < 1 in- 
dicates that coagulation is expected to set in above the 
critical temperature (i.e., W crossing the value 1 from 
above) within a rather narrow temperature interval. For 
the plots in Fig. [3] we have chosen parameter values which 
allow one to focus on the region of the onset of coagula- 
tion. As one can infer from the stability diagram (Fig. [2]) 
a = 10 corresponds to a relatively strong repulsion in 
the sense that for this value the range of values of E for 
which W < 1 holds and thus coagulation occurs is rather 
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narrow. This implies that one finds W < 1 only for in- 
termediate values of (. Such a case is represented by 
the red curve in Fig. [3] which corresponds to too = —2, 
a = 10, and s = 10. We find a clear increase of W for 
£ > 4.5; for larger £, IF becomes even larger than 1. 
The effect of varying too, a, or s relative to the red curve 
can be inferred from the brown, blue, and green curves, 
respectively Since the height of the potential barrier is 
proportional to s, stability is significantly enhanced for 
larger values of s (see the green curve (s = 50) in Fig. [3]). 
The value of £, for which In (W) turns from positive to 
negative (negative to positive) is larger (smaller) than 
the corresponding one for s = 10. For the composition 
corresponding to too = —5, which is closer to the crit- 
ical composition (see Eq. (|8|), the region of values of £ 
for which W < 1 is smaller than for too = — 2 (see the 
brown curve in Fig. [3j. Upon decreasing a, the value of 
£ for which In (W) changes from being positive to being 
negative becomes smaller (see the blue curve (a = 5)). 
For such small values of a, In (W) remains negative upon 
increasing £ (also for values of the parameters too and s 
not shown in Fig. 

In summary, along the thermodynamic paths consid- 
ered (and typically realized in experiments) the range 
of the attraction due to the CCFs grows steadily with 
increasing £, but the amplitude of the CCFs is a non- 
monotonous function of £. The CCF attains its maxi- 
mal strength for an intermediate value of £. Therefore, 
upon increasing £ the repulsive barrier in the effective 
pair potential is at first reduced; for sufficiently small 
values of a it can even disappear altogether. For even 
larger values of C, the potential barrier may grow again 
(or emerge) because the amplitude of the CCF decreases 
again. This latter barrier is located at smaller values 
of D than the one which is present for small values of £. 
The nonmonotonous dependence of the maximal strength 
of the CCF on temperature results in a nonmonotonous 
behavior of W = W (C)- Although, as stated above, a 
repulsive barrier may emerge again, it is likely that ag- 
gregates, which have formed for intermediate values of £, 
will not break up. This might be due to either specific 
microscopic interactions or due to the secondary mini- 
mum in the effective pair potential, which is still quite 
deep. Also in the case in which for all values of £ a re- 
pulsive barrier remains, coagulation can appear, due to 
a deep secondary minimum. This case will be discussed 
in Sect. EH 



V. STRUCTURE 

The bulk structure of a colloidal suspension is char- 
acterized by its radial distribution function g (r) which 
can be interpreted as the probability to find a colloidal 
particle a distance r apart from another particle fixed at 
the origin. The total correlation function (TCF) h (r) is 
related to the former one according to h (r) = g (r) — 1. 
The Ornstcin-Zernicke equation expresses h (r) in terms 



of the direct correlation function (DCF) c (r) and the 
number density p ) 51 i 52 

h(r)=c(r)+p J h (r') c(|r - r'|) dV, (9) 

or equivalently in Fourier space 

h(q)=c(q)/(l-pc(q)), (10) 

with h (q) = Je' qr /i(r) d d r, analogously for c(q), and 
q = |q|. In order to determine the two unknown functions 
c (r) and h (r) a second equation is required. This is the 
so-called closure equation which in its most general form 
is given by 

h (r) + 1 = cxp {-U (r) + h(r)-c(r) + b (r)} (11) 

where b (r) is the so-called bridge function, which in gen- 
eral is not known. The most common approximations 
for the bridge function are the so-called Percus-Yevick 
approximation (PY), 

b PY (r) = In [h (r) - c (r) + 1] - h (r) + c (r) , (12) 

and the so-called hypernetted-chain approximation 
(HNC), 

b H NC (r) = 0. (13) 

For practical reasons, in order to handle the hard core 
U (r < 2R) = oo, it is useful to introduce the function 

k(r) = h{r)-c{r). (14) 

In terms of k (r) the PY closure (Eq. (JT2J) ) can be written 
as 

cpy (r) = [exp {-U(r)} - 1] (fc (r) + 1) (15) 
and the HNC closure (Eq. ([13])) as 

chnc (r) = exp {-U (r) + k (r)} - k (r) - 1. (16) 

The correlation functions can be calculated 
iteratively i 52 ' 53 For a given approximate Ci (r) one 
calculates the TCF h l (r) according to Eq. (fTUf . By 
using a closure, i.e., by choosing a bridge function, 
Eq. (fTTj) renders a DCF c^+i which typically differs 
from Ci. This procedure is continued until satisfactory 
convergence is achieved. The initial guess c^ = i (r) is 
guided by the shape of the direct interaction potential. 

The applicability and reliability of this integral equa- 
tion approach (IEA) is discussed in detail in Ref. For 
comparison with Monte Carlo simulations in the case 
of a pair potential with attractive and repulsive parts 
see, e.g., Ref. Hll which discusses particles interacting 
with a pair potential containing attractive and repulsive 
Yukawa-like contributions e;e _Kir jr. The IEA is capable 
to reveal the rich phase behaviour of such systems. 
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1 2 3 4 5 6 
D/R 

FIG. 4. The radial distribution function g (r) as obtained 
from the HNC closure (Eq. JT3J) for four values of C ~ \t\~ v 
and for a volume fraction r\ = 0.04 of the colloids. The pa- 
rameters of the corresponding effective potential (attractive 
at large distances (la), Eq. l[5"|l) are a = 100 for the strength 
of the repulsion, mo = —100 for the reduced composition 
(Eq. ([8])), and s = k,R = 10; k _1 is the range of the soft 
repulsion, R is the radius of the colloids, and D = r — 2R is 
the surface-to-surface distance. The sequence of curves corre- 
sponds to the generic isochoral thermodynamic path realized 
in experiments. For these values of the parameters, upon de- 
creasing the reduced temperature t, i.e., increasing £, a single 
attractive minimum develops in the effective potential. This 
increasing attraction due to the CCF leads to the enhance- 
ment of g (r) close to the surface of the colloid. This peak 
suggests that the formation of colloid dimers is favored. The 
inset provides an enlarged view of the vicinity of the second 
maximum of g(r). 



The structure factor 

S(q) = l + ph(q) = 1/(1 -pc(q)) (17) 

can be determined by scattering experiments. From sum 
rules it follows that the isothermal compressibility xt is 
given byi£ 

WmS{q)= pk B T X T. (18) 

g-s-0 

The system becomes unstable for \t — > oo, correspond- 
ing to the critical point and, within MFT, to the spin- 
odals in the phase diagram (see, e.g., Ref. \53t) . 

The diversity of possible shapes of the effective pair 
potential given in Eq. © (see Fig. [T]) leads to a variety 




0.1 0.2 0.3 0.4 
D/R 

FIG. 5. The radial distribution function g (r) as obtained 
from the HNC closure (Eq. (|13[) ) for various values of £ ~ \t\~" 
and for a volume fraction r\ = 0.04 of the colloids. The pa- 
rameters of the corresponding effective potential (repulsive at 
large distances (lr), Eq. (0) are a — 40 for the strength of 
the repulsion, mo = —3 for the reduced composition (Eq. JHJl), 
and s = kR — 50; k," 1 is the range of the soft repulsion, R 
is the radius of the colloids, and D — r — 2R is the surface- 
to-surface distance. For these values of the parameters, upon 
decreasing the reduced temperature t, i.e., increasing £, an 
attractive minimum develops in the effective potential at in- 
termediate values of D while repulsion remains at small and 
large values of D. Due to this attraction a peak develops upon 
decreasing t. For a comparable potential well depth this peak 
is higher and narrower than in the case without a long-ranged 
repulsion (compare Fig. [4]). 

of radial distribution functions g (r) (see Figs. 2] and [5] in 
which the results for the HNC closure are shown; the re- 
sults obtained using the Percus-Yevick closure are almost 
the same, differing mostly in that, for the same value of 
C, the amplitude of the resulting g (r) is slightly smaller 
than the one within the HNC closure) . For temperatures 
far away from the critical temperature of the solvent (i.e., 
for ( = k£ <C 1), the colloids are effectively hard spheres 
with an effective diameter a > 2R due to the soft re- 
pulsive background contribution U rep . Accordingly, for 
such values of £, g (r) has the corresponding characteris- 
tics of a fluid of hard spheres, such as the rather broad 
first peak for small values of D. Due to the emerging 
attractive CCFs, for increasing ( ~ \t\ " , the radial dis- 
tribution function g(r) is enhanced close to the surfaces 
of the colloids. This implies an enhanced short-ranged 
order. 
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In the case that the effective potential exhibits a re- 
pulsive barrier at small values of D and is attractive 
throughout forge distances (denoted by (la) in the fol- 
lowing), we find that upon increasing £ the positions of 
the extrema of g (r) move towards smaller values of D 
(especially for intermediate values of Q and the value of 
g (r) at the maxima (minima) increases (decreases); the 
latter occurs mostly for larger values of C (see Fig. [2| . 

In contrast, in the case of effective potentials which 
are repulsive at forge distances (denoted by (lr) in the 
following), a new first peak at small values of D emerges 
and increases with increasing £ while the former first peak 
decreases and finally disappears (see the brown curve for 
C = 2.5 in Fig. [5] showing the emergence of a new peak 
at D/R ~ 0.1; the broad former first peak is located at 
larger distances D/R ~ 0.25). A similar trend is seen 
for the other extrema at larger D. However, since in 
that range of D the amplitude of g (r) is already quite 
small, the actual changes of the shape of g (r) cannot be 
described in such simple terms. 

In homogeneous, i.e., not phase separated systems, the 
(newly emerged or shifted) first peak in g (r) can become 
significantly larger in the (lr) case than in the (la) case 
(see Figs. [5] and [4] respectively, and note the different 
scales). The reason for this difference can be under- 
stood as follows. To a large extent the height of the 
first maximum in g (r) depends on the depth of the at- 
traction well of the effective potential. Indeed, for com- 
parable potential well depths the first peak in g (r) is 
only slightly larger in the (lr) case than in the (la) case. 
However, for equally deep potential wells the overall at- 
traction is stronger in the (la) case due to the attraction 
throughout large distances (which gives rise to a larger 
compressibility xt which is proportional to the integral 
of g (r) — 1). Accordingly, the systems interacting via 
(la) effective potentials can phase separate for shallower 
potential wells. Systems with potential well depths cor- 
responding to large first peaks in g(r) (as the ones shown 
for the (lr) potentials in Fig. O and with a (la) charac- 
ter would already be phase separated. In this sense, the 
repulsion occurring at larger values of D enhances the 
stability (w.r.t. phase separation) even of colloidal sus- 
pensions exhibiting large peaks in g(r) in homogeneous 
systems. In order to illustrate this point, we note that the 
compressibility corresponding to the blue curve (C = 2.6) 
in Fig. d] turns out to be approximately the same as the 
one corresponding to the green curve (( = 3.8) in Fig. [5l 
although for the latter system the peak in g(r) is nearly 
20 times larger. 



VI. COMPARISON WITH EXPERIMENTS 

For dilute suspensions, the second virial coefficient^ 
£?2 provides information about the strength of the radi- 
ally symmetric attraction between spherical particles: 
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FIG. 6. The value of the second virial coefficient B2 (Eq. (JT9J)) 
normalized by its value {B2) reg far away from Ti s ' as function 
of the temperature deviation from the critical temperature 
Tj s ' of the solvent and at its critical composition (E = 0). The 
same data in the inset but as function of the correlation length 

£ = £0 1 — T/Tc of the solvent. The experimental values 

(symbols) for latex particles immersed in the critical water- 
lutidine mixture are taken from Fig. 4 of Ref. [10(b). The 
lines correspond to the values obtained from the effective pair 
potential given in Eq. (JSJ. The values of a and s were obtained 
by using Eq. (|20l) for A and the values T = 5.7^C/cm 2 and 
2R = 555nm as given in Ref. 0(b). We have considered 
various values of k within the range reported in Ref. llClfb). 
We have adopted the values £0 = 0.25nm and T c = 310K. 
For the critical mass fraction ljl,c = 0.287 of lutidinei°' b ' one 
obtains 5 ^ e = 22.3. For further details see the main text. 



Beyond the ideal gas contribution it determines the lead- 
ing non-trivial term in the expansion of the pressure 
p (p) j (ksTp) = 1 + B2P + ... in terms of powers of the 
number density p. 

In Ref. [13 the second virial coefficient B2 has been de- 
termined for latex particles immersed in a mixture of 2,6- 
lutidine and water by using light scattering. In Fig.[6j for 
the critical composition (E = 0) we compare the varia- 
tion of Bi upon approaching the critical temperature T c 
for the effective pair potential in E q. O with the exper- 
imental data given in Fig. 4 of Ref. [13(b). In order to be 
able to compare these data with U (r) given in Eq. (fT|) we 
have to specify the amplitude A of the repulsive contri- 
bution U rep (r) = A exp {—tiD). This amplitude depends 
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on the surface charge density T of the colloids^ 

A=(ee y 1 T 2 K- 2 R/{k B T), (20) 

where e is the permittivity of the solvent relative to vac- 
uum, eo is the permittivity of the vacuum, and k is the 
inverse Debye screening length. We have adopted the ex- 
perimental values T = 5.7/iC/cto 2 and 2R = 555nm as 
given in Ref. Iiol(b) . The relative permittivity e as func- 
tion of the mass fraction uo L of lutidine can be obtained 
by using the Clausius-Mossotti relation^ b ^. We have 

(s) 

taken the values £o = 0.25nm and T c = 310K as typ- 
ical for water-lutidine mixtures (see, e.g., Ref. SI). In 
Ref. [Tol(b) no precise value of the Debye screening length 
is given; it is stated that « (7 ... 10) nm. There- 
fore in Fig. [6] we show results for B2 calculated for several 
values of ft; -1 [rim] = 5, 7, 10. 

Within their analysis, for temperatures far away from 
T c (s) (i.e., \T - T c (s) | > 3K) the authors of Ref. [13(b) 
obtained values for B2 which would correspond to a sys- 
tem of hard spheres of a radius of 4200nm. These ex- 
tremely large, unexpected values (they actually used par- 
ticles of diameter 555nm) cannot be easily explained by 
a soft, repulsive interaction^'-' 3 '' . In order to circumvent 
this non-trivial effect, in Fig. [6] we plot B2 normalized 
by the corresponding values {B2) reg far away from T c . 
For the normalization of the experimental data, we have 
used the mean values corresponding to the three largest 



T -T, 



which have been re- 



tcmpcraturc differences 
ported. For {B2) reg in the theoretical curves we have 
used the values as obtained by taking in Eq. ([5]) for x > 
only the soft, repulsive term into account (i.e., taking 
C = «£ = 0), so that (B 2 ) reg /v HS = 5.28,6.00,7.25 for 
/t _1 [nm] = 5,7, 10, respectively, where Vhs = (47r/3)i? 3 
is the volume of the hard spheres. 

Although we have adopted a rather simple form for the 
background repulsive potential, for k^ 1 = 5nm we obtain 
a fair agreement between the normalized experimental 
and theoretical values. Yet due to large error bars of the 
experimental data and large differences between not nor- 
malized experimental and theoretical values of B2 (the 
origin of which is unclear) we do not claim that this value 
of k^ 1 is the actual experimental one. Rather we consider 
in that case as an effective fitting parameter which 
accounts for differences between the experimentally re- 
alized and the model background potential in order to 
match the model with the experimental data. Further- 
more, we have not checked the sensitivity of our results to 
uncertainties in the other parameters. Both from experi- 
ment and theory it is evident that the attraction close to 
the critical temperature increases significantly, resulting 
in a steep decrease of B 2 towards large negative values. 

In Ref. fiolfb) the variation of B 2 as function of 

(s) 

T — Tc is reported also for a mass fraction uu L of lu- 
tidine larger than its critical value ojl,c — 0.287 (Fig. [7]). 
(This value of lo l , c = 0.287 is quoted in Ref. [13(b); it dif- 
fers slightly from the value uil, c = 0.286 quoted in Ref. S 
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FIG. 7. The variation of the second virial coefficient B2 
(Eq. ((T9J ) upon approaching the critical temperature Tc for 
an off-critical composition of the binary solvent. The same 
data in the inset, but as function of the correlation length 

£ = £0 1 — T/Tc \ of the solvent. The experimental val- 
ues (x) are taken from Fig. 2 of Ref. 10(b) and the three 
theoretical curves are obtained by using the effective poten- 
tial in Eq. (J5]). The values K~ 1 [nrn] — 7, 6, and 5 correspond 
to mo = -10.6, -13.5, and - 18.0, respectively (Eq. (©). 
For the lutidine mass fraction ujl = 0.35 > ujl,c the solvent 
is poor in the component preferred by the colloids. Therefore 
the critical Casimir force is stronger than for the critical com- 
position. Thus B2 can be expected to decreases more rapidly 
than for ujl = ujl,c- This is indeed the case for the theoretical 
curves. For the experimental data of Ref. 03(b) no depen- 
dence onui is seen; the dotted curves and the experimental 
data (o) correspond to the critical composition as shown in 
Fig. [6] B2 is normalized in the same way as in Fig. [fj The or- 
der parameter amplitude used here is — 0.765 (sec Eq. 
Ref. [73] in Ref.Q, and the text below Eq. (27) in Ref.0). By 
making use of the Clausius-Mossotti relation^ 13 we obtain 
e = 19.3 (see Eq. (|20[) ~). The values of the other parameters 
are the same as in Fig. [6] 



and as obtained from fitting the two-phase coexistence 
curve determined experimentally in Refs. and0 to the 
functional form \uo L — u>l,c\ = \t\^ with /3 = 0.3256; 
see Fig. [51 In other papers values are quoted within the 
range w L , c = 0.28... 0.29; see, e.g., Refs. II US, and SI) 
In general, for off-critical compositions the strength and 
the range of the critical Casimir force results from the 
interplay of the two bulk correlation lengths £ ~ \t\ v 
and £W ~ \h\ J \~ l '^^ S ^ (sec Eq. (|T|) and the discussion 
thereafter). For uo L = 0.35 the binary solvent is poor 
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in the component preferred by the kind of colloids used 
(i.e., water). In such a case the interplay of £ and £^ 
is rather complex due to the occurrence of capillary con- 
densation which leads to bridging between two colloids. 
In particular, at the same temperature the CCF can be- 
come stronger upon moving away from the critical com- 
position. But at off-critical compositions it is in general 
shorter ranged than at the critical composition. More- 
over, at off-critical compositions the CCF emerges as 
function of temperature only closer to the critical tem- 
perature as compared with the corresponding thermody- 
namic path at the critical composition. In Fig. [7] we com- 
pare the reduced experimental data for B2 corresponding 
to the mass fraction ui L = 0.35 > uil.c = 0.287 with the 
theoretical curves obtained within our effective model; 
for lo l = 0.35 one has {B 2 ) reg /vrs = 5.30, 5.65, 6.03 for 
k [nm] = 5,6,7, respectively. As can be inferred from 
Fig.[7]for = 5nm the total attraction (within our ef- 
fective model), in terms of B2, is stronger for uj l = 0.35 
than for the critical composition Wi iC , i.e., for the same 
value of £, B2 is more negative. For the theoretical curve 
corresponding to = 7nm the aforementioned inter- 
play of £ and manifests itself more clearly. For the 
off-critical composition the attraction due to the CCFs 
sets in only for larger values of £ but becomes more 
rapidly stronger with increasing £ so that, for large val- 
ues of £, B2 is more negative for lol = 0.35 than for lol,c- 
At the critical composition, B 2 starts to decrease (due to 
the attractive CCFs) already for smaller values of £ than 
for lo l = 0.35. 

For the experimental data (within the large error bars 
and the reported temperature range) no dependence of 



0.2 



on the com- 



thc value of B2 as function of T — T< 
position is observed. (Reference llOKb) provides the val 



ues of B2 as function of 



T - 



but lacks an absolute 



■.(»). 

L CX 1 



t£> is the 



value for the coexistence temperature T c 
temperature at which the solvent undergoes a first-order 
phase segregation for to L = 0.35 ^ Wi jC ; see Fig. 1(a) in 

Ref. 0. We take T$ - T c (s) = 0A45K as obtained^ from 
|wl -u L ,c\ = &u \tf with B u = 0.765.) This finding 
is in contrast to the behavior of the measured effective 
particle- wall interaction potentials of Ref. [13, for which 
a dependence on the composition of the water-lutidine 
mixture has been observed. At present it remains un- 
clear whether the dependence of B2 onwi could not be 
resolved experimentally or whether the data from these 
experiments happen to lie in that crossover regime in 

which the curves B2 — T^j happen to be close to 

each other for the critical and the off-critical mixture 
(see the inset in Fig. [7| . 

For the same system the authors of Ref. fiol(b) used 
for the measurement of B 2 , in Ref. 0(b) the aggregation 
phenomena of immersed colloids have been studied. We 
have used the corresponding given experimental values of 
the parameters, and especially Eq. ([2"U|) for the repulsion 
strength parameter A, in order to investigate these state 



phase segregated 
region 




FIG. 8. Comparison of the values of B2 (Eq. ((T9|) ') at the 
thermodynamic state points of the solvent for which immersed 
colloids aggregate. The four experimentally obtained state 
points (squares; the straight black dotted lines in between 
are a guide to the eye) are taken from Fig. 1 (middle) of 
Ref. |7|(b)). The isolines of constant B2 (full colored lines), 
each belonging to one of the state points (squares), are cal- 
culated within our model (Eq. l[5]l). Each i?2-isoline can cap- 
ture some qualitative trends of the possible shape of the line 
of the onset of aggregation, but they cannot serve as quan- 
titative indicators for the onset of aggregation. For further 
details see the main text. In order to obtain the theoreti- 
cal values of B2 we have taken £0 = 0.25nm, = 6nm, 

(s) 

Tc = 307K, and lol, c = 0.286 and we have used the exper- 
imentally given values 2R = 378nm and T = 3.85(iC/cm 2 . 
The repulsion strength A is obtained from Eq. (|20[) with e 
depending, according to the Clausius-Mossotti relatior*22 b *££, 
on the mass fraction lol of lutidine. The experimentally de- 
termined coexistence points (circles) on the binodal of phase 
segregation with a lower critical point agree well with the rela- 
tion \lu l - lu l , c \ = B u \tf, where j3 = 0.3256 and B w = 0.765 
(gray line); in accordance with a lower critical point the re- 
duced temperature is t 



(rj a 



T) IT, 



points within our model. 

First, it is worthwhile to mention that for the param- 
eters corresponding to these experiments, the repulsive 
barrier for small values of x = kD is of the order of 
tens of ksT. Therefore, it seems that the observed ag- 
gregation is due to the very deep secondary minimum 
emerging at a certain distance D m i n . This is a differ- 
ent mechanism than the one discussed in Sect. HVl where 
the aggregation due to the disappearance of the repul- 
sive barrier at short distances D has been considered. 
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While in the latter case the colloids would stick together 
with their surfaces at contact, in the former case they 
are close but not at contact. With K,D m i n » 1 and 
s = kR w 10 . . . 100, from the effective potential given 
in Eq. ([5]) it follows that the colloids are a surface-to- 
surface distance D/R w 10 _1 . . . 10" 2 apart. In Ref. 0(b) 
the mean surface-to-surface distance has not been esti- 
mated, but in experiments performed for a similar sys- 
tem (see Ref. @) it has been found that in the floccu- 
lated and sedimented aggregates the average interparti- 
cle distance is about three times the mean radius of the 
polydisperse colloids. This is in the stark contrast with 
our theoretical estimates. Other experimental reports 
do not specify the average interparticle distance occur- 
ring in the aggregates. In Ref. [lj it is stated, that the 
sediment is less dense than expected for a tight pack- 
ing of spheres and in Ref. Il7| the aggregates are char- 
acterized as being compact, not fractal. The authors of 
Ref. [TH argue that for their system the aggregates are in- 
deed crystalline. The possibility, that the flocculated col- 
loidal particles form a bona fide, thcrmodynamically sta- 
ble, liquid-like or crystalline phase cannot be addressed 
within the effective, one-component approach used here, 
because within this approach the reliable prediction of 
the occurrence of a phase transition can be made only 
if the colloidal number density is sufficiently large (see 
Ref. 0, especially, Subsect. Ill B 3), whereas in the ex- 
periments referred to above the packing fraction is as 
small asij w 10" 6 ...10- 3 . 

Next, we address the issue whether a relation can be 
established which connects the onset of the reversible ag- 
gregation to a threshold of a specific quantity which is 
both theoretically and experimentally accessible. To do 
so we follow Ref. [l(J and analyze whether the values of B 2 
corresponding to thermodynamic state points of the sol- 
vent, at which the aggregation of colloids sets in, are sim- 
ilar. For such state points, as determined experimentally 
and reported in Fig. 1 (middle) of Ref. 0(b), we have cal- 
culated the corresponding values of B 2 within our model 
for four values of uj l yielding^ B 2 / (^-R 3 ) = -7400 
(red square), —560 (blue), —200 (green), and —76.0 (yel- 
low) (Fig. [5]). In Fig. [5] the isolines of constant B 2 for 
these values of B 2 are shown (using the same color code). 
Clearly the four curves do not coincide to form a single 
curve with the four squares on it which would support 
the hypothesis that aggregation sets in at a threshold of 
value of B 2 . However, one has to keep in mind that the 
CCFs strongly depend on the temperature. This depen- 
dence is amplified within the resulting value of B 2 due 
to its exponential dependence on the potential. There- 
fore the value of B 2 may vary significantly even within 
a narrow temperature interval (compare Figs. [6]and[7|. 
Moreover, small uncertainties in the measured tempera- 
ture, at which the aggregation sets in, result in a large 
variation of B 2 . Each £?2-isoline reproduces some of the 
qualitative trends of the experimental data points for ag- 
gregation onset, but there is no satisfactory quantitative 
agreement. Either the isoline is too narrow (see the red 



curve which cannot capture any of the yellow, blue, or 
green squares) or the isoline extends to too large values 



of 



T-T, 



0) 



(see, e.g., the yellow curve which misses 

the blue, red, and green squares). Besides the possibil- 
ity that B 2 is indeed not the same at the various state 
points at which coagulation sets in, there are further rea- 
sons which could explain this disagreement. One reason 
could be the simplified choice for the background poten- 
tial which we have used. Furthermore, there are several 
experimental parameters for which we could use only es- 
timates. Within their analysis, the authors of Ref. 0(a) 
conclude that for calculating the electrostatic potentials 
one should not use the bare surface charge density as 
given by the manufacturer. Rather, the dissociation of 
the colloidal surface in the solvent has to be taken into 
account. In our analysis we took the nominal value for 
the surface charge due to the lack of more detailed data. 
The authors of Ref. point this out, too, and reveal the 
role of impurities as another important factor. They de- 
termined the aggregation line in "fresh" mixtures and in 
the same mixture but one or several days later. While 
the aggregation phenomenon as such was reproducible, 
the values of the thermodynamic variables of the solvent 
at which the aggregation sets in varied strongly within 
one and the same sample. 

In the same reference, a dependence of the aggrega- 
tion line on the volume fraction 77 of colloids is observed, 
whereas in Ref. UH it is stated that the shape of the region 
in thermodynamic space where aggregates are forming 
does neither depend on the colloidal concentration nor 
on the colloidal radius. In Ref. 0(b) a dependence of the 
time scale of the aggregation process on rj is reported. 
Our analysis does not capture such a dependence. 

In conclusion, our results suggest that the critical 
Casimir force plays an important role in order to ex- 
plain the aggregation phenomenon in this region of the 
phase diagram of the solvent. Yet, in order to be able 
to determine all relevant mechanisms and forces and in 
order to obtain quantitative agreement between the ex- 
periments and the theoretical results, for both further 
efforts are needed to provide more precise values of the 
relevant quantities. 

Finally, in Fig. [5] experimental data for the binodal of 
the water-lutidine mixture are shown, which nicely agree 
with the functional form ui^ —ujl,c = ±£L \t\ (gray line), 
where p = 0.3256 and = 6.765. This shows that 
our description, which is symmetric w.r.t. ujl = ujl,c 7 is 
sufficient for the phase space region of interest here. 



VII. SUMMARY 

We have studied aggregation phenomena of monodis- 
perse colloidal suspensions with near-critical solvents. 
In such solvents, driven by the diverging bulk correla- 



te 



0= 



tion length t(t= (V c (s) - T) /T< 
so-called critical Casimir forces act on the colloidal par- 
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tides, in addition to regular background forces. In order 
to obtain our results, we have adopted a description in 
terms of an effective one-component system. Within this 
approach the presence of the solvent is taken into account 
via contributions to the effective pair potential of the col- 
loids. We have considered an effective pair potential for 
the colloids (with radius R) which consists of the critical 
Casimir contribution plus a regular background poten- 
tial. For the latter we have employed a generic, softly 
repulsive contribution with strength A and acting on a 
length scale We have focused on dilute suspensions 
with low colloidal number densities p, for which the ap- 
proximation of pairwise additive critical Casimir forces 
(CCFs) is valid. In addition, we have applied the Der- 
jaguin approximation (Eq. (|2|l). within which the pair 
potential between two spheres can be expressed in terms 
of the scaling function of the CCF in a slab. In order to 
obtain access to the dependence of the universal scaling 
function for the CCF on the scaling variable £ associ- 
ated with the bulk ordering field hi, in spatial dimension 
d = 3 we have employed an approximation based on a 
combination of both Monte Carlo simulation data and of 
mean-field theory results (Eq. (J3])). The resulting effec- 
tive pair potential is given in Eq. © and various exam- 
ples are plotted in Fig. [TJ On this basis we have obtained 
the following main results: 

1. We have calculated the stability ratio W ~ 
tpair/tdiff (Eqs. ((6]) and ((7|)), which is a measure 
of the stability of a colloidal suspension. The time 
scales tp a i r and tdiff are related to the pairing of 
two colloids and to the diffusion of a single particle, 
respectively. Rapid coagulation is expected to set 
in for W < 1. In Fig. [5] the contour lines W = 1 
are shown for various values of £ = k£ ~ |i| 
and s = kR in the space spanned by the parameter 
a = A/s and by the scaling variable S. For the case 
of a binary liquid mixture £ is related to the differ- 
ence of the chemical potentials of the two species. 
This shows that the CCF can lead to rapid coag- 
ulation. The coagulation region extends to larger 
values of a for £ < because the CCFs are stronger 
for compositions of the solvent slightly poor in the 
component preferred by the colloids. 

2. The variation of W with temperature at a fixed 
solvent composition, as it is typically realized in 
experiments, is shown in Fig. [3] We find that W 
rapidly decreases to W < 1 upon approaching the 
critical temperature via the disordered phase. This 
implies, that the aggregation should set in within 
a narrow temperature interval. 

3. Concerning the structure of the colloidal suspen- 
sion, we have solved numerically the Ornstein- 
Zernicke equation (Eq. ([9])) using the Percus-Yevick 
and the hypernetted-chain closure (Eqs. (|T5j) and 
(HHJ), respectively). In Figs. 2] and [5] we show the 
radial distribution function g (r) for two different 



parameter sets (a,s,mo) and various values of £, 

where m^ 1 ~ c a — Ca,c is related to the composi- 
tion of the solvent (Eq. ©). Close to the surface 
of a colloidal particle (r > 2R) a peak develops 
for increasing values of reflecting the increasing 
attraction due to the CCF. This indicates the pres- 
ence of a pronounced short-ranged order and that 
the formation of colloidal dimers is favored. 



4. In Figs. [6] and [3 for the critical and for an 
off-critical composition, respectively, we compare 
the variation of the second virial coefficient B2 
(Eq. (HHJ)) upon increasing £ = £0 \t\~", as obtained 
within our model, with experimental data^^ . For 
both, the theoretical results and the experimental 
data, for increasing £ the virial coefficient B2 de- 
creases towards large negative values. Due to ex- 
tremely large experimental values of I-B2I also far 

(s) 

away from T c , the origin of which is not clear, the 
quality of the quantitative comparison is limited. 
For the off-critical composition, the CCF can be- 
come stronger than at the critical composition but 
it has in general a shorter range and occurs closer 
to the critical temperature of the solvent. There- 
fore, for the same value of £, the total attraction 
due to the CCFs can either be smaller or larger 
for an off-critical composition compared with the 
critical one. Accordingly, as can be seen in Fig. 
B2 can be less or more negative, respectively. For 
the experimental data of Ref. Hol(b) there is no vis- 
ible dependence on the composition of the solvent, 
which could be due to this crossover behavior. 



5. We have addressed the issue whether a relationship 
can be established between the onset of aggregation 
and a specific quantity which is accessible both the- 
oretically and experimentally. To this end, within 
our model we have calculated B 2 at state points 
for which in experiments aggregation sets in. The 
corresponding £>2-isolincs capture qualitatively the 
main trends of the aggregation line, but no satis- 
factory quantitative agreement has been achieved 
(see Fig. [8]). We have discussed various possible 
explanations for this discrepancy. 



To conclude, our results show that CCFs, which can be 
easily controlled by temperature T and the strength and 
range of which can be varied by changing the bulk order- 
ing field hi,, can lead to aggregation in colloidal suspen- 
sions with near-critical solvents. Within the approach of 
using an effective potential, we have identified the ranges 
of values for the background repulsive potential and the 
values of the scaling variables associated with the critical 
solvent, for which the colloids aggregate. 
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parameters in the potential U (r). For the cases considered here a j (2R) varies within the range 1.03 . . . 1.05. 



